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Interacting two-component Fermi gases loaded in a one-dimensional (ID) lattice and 
subjected to an harmonic trapping potential exhibit interesting compound phases in 
which fluid regions coexist with local Mott-insulator and/or band-insulator regions. Mo- 
tivated by experiments on cold atoms inside disordered optical lattices, we present a 
theoretical study of the effects of a correlated random potential on these ground-state 
phases. We employ a lattice version of density-functional theory within the local-density 
■ approximation to determine the density distribution of fermions in these phases. The 

f ' exchange-correlation potential is obtained from the Lieb-Wu exact solution of Fermi- 

' Hubbard model. On-site disorder (with and without Gaussian correlations) and harmonic 

f^*^ , trap are treated as external potentials. We find that disorder has two main effects: (i) it 

destroys the local insulating regions if it is sufficiently strong compared with the on-site 
atom-atom repulsion, and (ii) it induces an anomaly in the inverse compressibility at low 
KJ{ density from quenching of percolation. For sufficiently large disorder correlation length 

I ^N , the enhancement in the inverse compressibility diminishes. 

Keywords: Fermi-Hubbard model; optical lattices; disorder. 



1. Introduction 

Disorder and interaction effects in condensed matter systems have a long and rich 
history. The interplay between them has been the subject of continuing interest. The 
notable examples range from two-dimensional electron systems with long-ranged 
Coulomb interactions^ to liquid 4 He absorbed in various substances such as aerogel 
and vycor^and granular superconductors.^ In the former system experimental and 
theoretical investigations^ reveal the crucial role played by disorder as the two- 
dimensional electron system undergoes a metal-insulator transition. Furthermore, 
thermodynamic quantities carry the signature of phase transitional^] along with 
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transport properties. 

In recent years, cold atomic systems are being used to investigate the interplay 
between single-particle randomness and interaction effects with the help of optical 
lattices allow ing access to strong coupling regimes through the depression of kinetic 
energy.^EE] Furthermore, the on-site interaction can be tuned either indirectly by 
changing the strength of the lasers that create the optical lattice or directly by 

iini 

means of a Feshbach resonance!^ 

After the initial observation of superfluid to Mott insulator trans ition using 
bosonic atoms in a n optic al lattici^ and creation of disorder potentials^] studies 
on fermionic atoms^H^ are beginning to be explored. A couple of forthcoming 
reviewi^l encompass many aspects of the theory of ultracold Fermi gases. 

In this work we study the interaction and disorder effects on a one-dimensional, 
two-component Fermi gas trapped in a harmonic confinement potential and an 
optical lattice. We use the lattice version of density functional theorjEl taking 
advantage of the exact solution of the one-dimensional Fermi-Hubbard model. 
Ground-state calculations in the absenc e of disorder have been performed within 
a number of numerical technique ^ 1 18 | to identify various phases. In the present 
approach, the harmonic confinement potential and the disorder potential are treated 
as part of the Kohn-Sham potential within a local-density approximation which has 
an exact representation for the uniform case. 

Generalizing our earlier worlff^ on the effects of uncorrelated disorder, we here 
study correlated disorder which depends on the correlation length as a parameter. 
The site occupation profiles and thermodynamic stiffness are determined to study 
the effects of disorder. Our motivation for considering correlated disorder comes 
from the experiments^ on Bose-Einstein condensates in optical speckle potentials 
where it was noted that the correlation length is several times longer than the lattice 
spacing. 

The paper is organized as follows. In Sec. 2 we briefly outline the Fermi-Hubbard 
model in a harmonic potential and disorder. We then give the main ingredients of 
our density-functional theory approach. Section 3 presents our results on site occu- 
pations in the presence of correlated disorder. We conclude with a brief summary 
in Sec. 4. 



2. Theory and Model 

We consider a two-component Fermi gas with N atoms constrained to move under 
harmonic confinement of strength V% inside a disordered ID optical lattice with 
unit lattice constant and L lattice sites i G [l,L]. The system is described by a 
single-band Hubbard Hamiltonian, 

L-l L 

K = -t XMo-Ci+l* + H.c.) + U^2 n lT n 4 

2—1 (J i—1 
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L L 

+ V 2 Y / ^-L/2fn t + Y / ^- (!) 
i=l i=l 

Here iy = t > if «, j are nearest sites and zero otherwise, a =|, J. is a pseudospin- 
1/2 label for two internal hyperfine states, r£i CT = cJ CT Cio- is the pseudospin-resolved 
site occupation operator, and hi = ^2 a h ia . The effect of disorder is simulated by 
the last term in Eq. ([1]). The above Hamiltonian without the confining potential 
and disorder term has been solved exactly using the Bethe Ansatz technique by 
Lieb and Wuff^ The exact solution provides us with the ground-state energy at 
any coupling strength u = U/t and filling n = N/L. In the noninteracting and 
unconfined limit (i.e., without the U and trap terms) one recovers the Anderson 
models studied intensively for the Anderson localization problem. The successful 
creation of optical lattices to study cold atomic systems has led to renewed interest 
in these model systems, since the atomic systems offer a wide range of control on 
the various parameters. 

In this work we generalize our previous study^^ on uncorrclated noise to Gaus- 
sian correlated disorder, defined as 



(2) 



2i 



where £ is the correlation length (in units of the lattice spacing, which is set to 
unity throughout this work) and Wj is randomly chosen at each site j from a uni- 
form distribution in the range [— W/2, W/2]. Because of the following mathematical 
identity, 

: ex P ^ — = Oij , (3) 

Eq. @ tends smoothly to uncorrelated (i.e. white) noise in the limit £ — > 0. We 
illustrate typical behavior of £j over a lattice of L = 200 sites in Fig. 1 for some 
values of £. As the correlation length increases the correlated disorder potential 
becomes smoother with amplitude smaller than that in the uncorrelated case. 

A particular set of values e% is a realization of disorder. Each realization defines 
an external potential V, = V%(i — N s /2) 2 + E{ which is the sum of a harmonic trap 
potential and disorder, the site-occupation functional theorjffSl (SOFT) to deter- 
mine the site occupation n, = (^|ni|$) where $ is the ground-state of Tt for this 
particular disorder realization. The site occupation A/i is obtained by the disorder 
ensemble average, A/i = {{rii))dis- SOFT is the discrete or lattice version of the den- 
sity functional theory which has been successfully applied to the present system in 
the absence^Sl and presencd^ of uncorrelated disorder. In the clean limit the local- 
density approximation which we adopt has been shown to be reliable through exten- 
sive comparisons with accurate quantum Monte Carlo calculations Local-density 
approximation based density-functional schemes for disordered systems have been 
employed to study the low-density compressibility anomaly in the two-dimensional 
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i 

Fig. 1. (color online) The correlated disorder potential e, as a function of the site position for a 
lattice of L = 200 sites. The disorder strength is W/t = 3. Different values of £ are indicated in 
the legend. 



metal-insulator transition (see Ref. [51 [5T]) and the statistical properties of 2D dis- 
ordered quantum dotsP^l 

The total energy is a unique functional of n,; which can be written as 

f[n]=5^e(n,u) + 2V i (2 i )n(z i ), (4) 

i i 

where e(n, u) is the ground-state energy of the Hubbard Hamiltonian as obtained by 
Lieb and Wuff^ The Euler-Lagrange equation that follows from the above energy 
density is 

de 

-g^ + VKs(zi) = constant , (5) 
where the Kohn-Sham potential is 

VKS = ^Un + v xc (z i )+Vi(zi), (6) 
Within the local-density approximation, exchange-correlation potential is approxi- 
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mated by 

R r 1 J 

(7) 



e(n, u) — e(n, 0) — ~^Un 2 



The above set of equations allow us to obtain the site occupations rij and their 
disorder averages Mi ■ 



3. Results and Discussion 



W= 
W= 1 




180 



i 

Fig. 2. (color online) Site occupation Mi as a function of site position i for N = 60 fermions 
with u = 4, £ = 4, and V^/t = 2.5 X 10 -3 in a lattice with L = 200 sites. The four curves have 
been calculated for different values of disorder strength: W/t = (solid line), W/t = 1 (circles), 
W/t = 4 (triangles), and W/t = 20 (squares). 



Before we describe the effects of disorder on the trapped lattice fermions, we 
outline the variou s gr ound-state phases of a clean system as obtained by previous 
numerical studiesff^ There are altogether five phases (A...S) controlled by the 
interaction strength U /t, number of fermions N , and number of lattice sites L. 
Phase A is a fluid with < n, < 2. In phase B a Mott insulated occupies the central 
region of the trap with rij = 1. In phase C a fluid with 1 < < 2 is embedded in the 
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Mott plateau. Phase I? is a band insulator with rij = 2 surrounded by fluid edges 
and embedded in the Mott plateau. Finally, in phase E a band insulator in the 
central region of the trap coexists with fluid edges. The sketch of site occupations 
in these phases was given in Ref. 19. 



W= 
W= 1 




180 



Fig. 3. (color online) Site occupation Mi as a function of site position i for TV = 70 fermions 
with u = 4, £ = 4, and V%/t = 2.5 X 10 -3 in a lattice with L = 200 sites. The four curves have 
been calculated for different values of disorder strength: W/t = (solid line), W/t = 1 (circles), 
W/t = 4 (triangles), and W/t = 20 (squares). 



In Fig. 2 we show the disorder-averaged site occupation Mi for a system of 
fermions with N = 60 in an optical lattice with L = 200 sites. The interaction 
strength is u = 4 and the trap potential is V^/t = 1.5 x 10 -3 . Here we fixed the 
disorder correlation length to be £ = 4 and varied the disorder strength W/t. The 
clean system is in phase B. As in the case of uncorrelated disorder-^ we find that 
with increasing W/t the Mott insulating region is depleted and the site occupation 
profile Mi broadens. However, the rate of depletion and broadening are smaller than 
the uncorrelated disorder. In other words, the Mott insulating region is more sta- 
ble against the formation of a disordered fluid phase in the presence of correlated 
disorder. 
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Fig. 4. (color online) Site occupation Mi as a function of site position i for N = 200 fermions 
with u = 8, £ = 4, and Vi/t = 2.5 X 10~ 3 in a lattice with L = 200 sites. The five curves have been 
calculated for different values of disorder strength: W/t = (solid line), W/t = 3 (solid circles), 
W/t = 5 (triangles), W/t = 10 (squares), and W/t = 20 (empty circles). 



In Fig. 3 we display the disorder averaged site occupation Mi for a system with 
N = 70 atoms keeping the rest of the parameters the same as in Fig. 2. N = 70 
is the critical number of atoms at which the phase transition B — > C occurs in the 
clean limit. For weak uncorrelated disorder a fluid phase with Mi > 1 is induced 
at the center of the trap. Essentially the same behavior is observed for correlated 
disorder (£ = 4). Sufficiently strong W/t destroys the Mott plateau completely 
and the system becomes a disordered fluid with Mi < 1. We find that the critical 
disorder strength for this to happen is much larger for correlated disorder than the 
corresponding value of W/t for uncorrelated disorder. 

In Fig. 4 we show the site occupation for a strongly interacting system (U/t = 8) 
with correlated disorder (£ = 4). The disorder free system is in the T> phase discussed 
above. With increasing W/t the band insulating region at the center of the trap is 
depleted. At the same time, the Mott insulating regions are also slowly destroyed. 
To see more clearly the rate at which the band and Mott insulating regions are 
destroyed, we plot in Fig. 5 the number of consecutive sites iVjviott and -/Ve an d such 



I 



I 



May 27, 2008 0:3 WSPC/INSTRUCTION FILE ldcorre 



8 Gao et al. 




Fig. 5. (color online) The number of consecutive sites iVMott and A^Band such that \Afi — 1| < 
and \Mi — 2| < 0, respectively, as a function of W/t. The parameters are the same as in Fig. 4. 



that \J\fi — 1| < and \J\fi — 2| < 0, respectively, as a function of W/t. Comparison 
with the same calculation in the uncorrclated disorder cas reveals the increased 
stability range of insulating regions for correlated disorder. 

The effect of disorder on ID fermions and in particular transitions between 
different phases can also be assessed through the thermodynamic properties. For this 
purpose we calculate the inverse compressibility defined as k" 1 = ((N 2 5[j,/5N))ai s 
where fi is the chemical potential. In Fig. 6 we show function of N for 

various values of the disorder correlation length £. The phase transitions can be 
identified as sharp kinks in a clean system (W/t = 0). Uncorrelated disorder has two 
main effects on k^ 1 as discussed in our previous workff21 Firstly, the sharp features 
indicating phase transitions are smoothed out. Secondly, a large enhancement of 
k^ 1 at low density is observed. This is reminiscent of a similar behavior found in 
2D electron systems. As N decreases the atoms mostly occupy the deepest valleys in 
the disorder landscape, thus the high density regions in the system tend to become 
disconnected. For a given interaction strength u and low N the system stiffens as 
disorder grows. In the present situation we observe that with increasing correlation 
length £ the low density enhancement in k _1 is diminished. The basic reason for 
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Fig. 6. (color online) The inverse compressibility k,~ 1 (in units of t) in the presence of correlated 
disorder as a function of N for V 2 /t = 2.5 X 10~ 3 , W/t = 5, u = 8, and L = 200 lattice sites. The 
correlation length £ = (circles), £ = 6 (squares), and £ = 32 (triangles). 

this behavior is that increasing £ makes the disorder potential more smooth and 
more shallow compared to the £ = case. Therefore the atoms are less localized in 
a correlated disorder potential. 

4. Concluding remarks 

We have studied the one-dimensional Fermi-Hubbard model in a harmonic con- 
finement potential and in the presence of correlated disorder. This is believed to 
represent cold fermionic atoms in an optical lattice created by standing laser waves. 
Our numerical calculations of the site occupations are based on a lattice version 
of density functional theory in which we make use of the exact solution of the 
one-dimensional Hubbard model to treat the exchange-correlation effects. Disorder 
affects the ground-state phases of the interacting Fermi gases confined in a har- 
monic potential and an optical lattice. The insulating regions appear to be stable 
against both uncorrelated and correlated disorder. The anomalous enhancement of 
the stiffness observed at low density for uncorrelated disorder decreases with increas- 
ing correlation length. We hope that our results on disorder effects will stimulate 
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experimental investigations with cold Fermi atoms in the future. 
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